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Abstract 

Using the isospin dependent quantum molecular dynamics model, we study the effect of charge 
asymmetry and isospin dependent cross-section on d((Cos2(f>)) ptdpt ' Simulations have been 
carried out for the reactions of X m + 124 X m , where m = (47, 50 and 59) and 40 Si6 + 40 S±q. 
Our study shows that these parameters depend strongly on the isospin of cross-section and charge 
asymmetry. The distribution of nucleons and fragments is not symmetric around the beam axis. 
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I. INTRODUCTION 



It is well known that collective flow is an important observable in heavy ion collisions 
(HIC) and it can give some essential information about the nuclear matter, such as the 
nuclear equation of state l|-|4|. For the last few years, collective flowhas been used as a 
powerful tool to explore the nuclear equation of state (EOS) as well as in medium nucleon- 



nucleon cross-section 



5] Anisotropic flow is defined as the different n th harmonic coefficient 



v n of the Fourier expansion for the particle invariant azimuthal distribution 



6]: 



A N 

— = l + 2Y,v n Cos(n<p) (1) 

n=l 

where 4> is the azimuthal angle between the transverse momentum of the particle and the 
reaction plane. Note that the z-axis is defined as the direction along the beam and the 
impact parameter axis is labelled as x-axis. Anisotropic flows generally depend on both 
particle transverse momentum and rapidity, and for a given rapidity the anisotropic flows 



at transverse momentum p t (pt=\/(pl +P«)> where p x and p y are projections of particle 



transverse momentum in and perpendicular to the reaction plane, respectively. The first 
harmonic coefficient V\ is called directed flow parameter. Directed flow is the measure of 
the collective motion of the particles in the reaction plane. This flow is reported to diminish 
at higher incident energies due to the large beam rapidity 7|. The second harmonic 



coefficient Vi is called the elliptic flow parameter u 2 . Elliptic flow in heavy ion collisions 



is a measure of the azimuthal angular anisotropy of particle distribution in momentum 
space with respect to the reaction plane 8j . The elliptic flow at intermediate energy HIC is 
complex phenomenon because it is determined by the interplay among fireball expansion, 
collective rotation, the shadowing of spectators, Coulomb repulsion, and so on. Both 
the mean field and two-body collision parts play important roles: the mean field plays a 
dominant role at low energies, and then gradually the two-body collisions become dominant 
with energy increase. The transverse radial dependent transverse velocity can reflect 
the correlation between spacial and momentum coordinates, and reveal the force change 
on fragments along the transverse radius. The magnitude of the elliptic flow depends 
on both initial spatial asymmetry in non-central collisions and the subsequent collective 
interactions. Experimentally observed out-of-plane emission termed as squeeze-out was 
observed by SATURNE (France) by DIOGENE collaboration |g|. The Plastic-Ball group 
at the BEVELAC in Berkley were the first one to quantify the squeeze-out in symmetric 
systems llOfl. The elliptic flow is sensitive to the properties of the dense matter formed 

n fl 

during the initial stage of heavy ion collision [llj and parton dynamics [12f at Relativistic 
Heavy Ion Collider (RHIC) energies. 



C. Pinkenburg et al., [l3j have measured an elliptic flow excitation function for midcentral 
collisions of Au + Au at 2, 4, 6, and 8 GeV/nucleon respectively. The excitation function 
exhibits a transition from negative to positive elliptic flow with transition energy E trans = 
4A GeV. 



J. Lukasik et al., 



2| studied the distributions of the squeeze angle for incident energies from 



3 



40 to 150 MeV/nucleon, their study revealed that, the minima at <fi = tt/2, occur at lower 
energies and more peripheral impact parameters while peaks at (p — 7r /2, most strongly 
pronounced in the more central bins at the higher incident energies. 

J. H. Chen et al, [4] studied the azimuthal angular distribution of raw <fi yields with 
respect to the event plane, they showed that, the finite resolution in the approximation of 
event plane as reaction plane smears out the azimuthal angular distribution and leads to a 
lower value in the apparent anisotropy parameters. In this paper our aim is to check that, 
whether the distribution of nucleon contributing to elliptical flow are distributed equally 



along the ellipse or not. Our study is performed within the framework of IQMD 14J model 



which is the improved version of QMD 



15] model. 



II. RESULTS AND DISCUSSION 



For the present analysis, simulations are carried out for the reactions of X m + X m , 



where m = (47, 50 and 59) and Sie+ <Si6- The phase space generated by the IQMD model 



has been analyzed using the minimum spanning tree (MST) 



16] method. The elliptical flow 



is defined as the average difference between the square of x andy components of the particles 



transverse momentum. Mathematically, it can be written as 



5]: 



(v 2 ) =< Cos2<p >= (4t4> 



(2) 



Px + Py 

where p x and p y are the x and y components of the momentum. The positive value of 
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FIG. 1: (color online) Azimuthal angle dependence of wr^^Wj > ^ or ^ ree nuc l eons (upper panel), 
LMF's (middle) and IMF's (lower panel). 

elliptical flow describes the eccentricity of an ellipse-like distribution and indicates in-plane 
enhancement of the particle emission. On the other hand, a negative value of v 2 shows the 
squeeze-out effects perpendicular to the reaction plane. Obviously, zero value corresponds 
to an isotropic distribution. 

To study the effect of isospin dependent cross-section and charge asymmetry on d ^Q^ 2 <j>)) > 
we display in fig. 1, the azimuthal angle dependence of d ^^ 2 ^)) > ^ or ^ ree nuc l eons (A = 1) 
(upper panel), LMF's (2 < A < 4) (middle) and IMF's (5 < A < A tot /Q) (lower panel) at 
an incident energy E = 100 MeV/nucleon for the reactions of 124 Ag^ 7 + 124 Ag i7 (left panels) 
and 124 Pr 59 + 124 Pr 59 (right panels). Figure reveal: 

(a) Minima at 2<fi=7r/2 indicate predominantly in-plane emission, while the peaks at 2<p= 



5 



and 7r, corresponds to a preference for azimuthal emission in-plane and perpendicular to 
the reaction plane, the so-called squeeze-out. 

(b) Peak is more pronounced at 2(p= than at 2<p= n, which indicates that number of 
particles emmited in-plane are large as compare to the number of particles emitted out-of- 
plane. This means that ellipse formed is not symmetric around the Z-axis i.e in the collision 
of symmetric nuclei the distribution of nucleon after the collision in momentum space is not 
uniform. 

(c) There is a very little influence of charge asymmetry on the variation of d ^^ 2 4>)) wr th 
(Cos2(f)}. d ({g o ^ 2 0)) i s more for neutron rich system 124 Ag 47 + 124 Ag 47 than neutron deficient 
system 124 Pr 59 + 124 Fr 59 due to increase in repulsive forces. 

(d) d{(Cos2<j>)) * s sens itive to different nucleon-nucleon cross-sections. Its value is more in case 
of isospin dependent cross-section. This happens because in the case of isospin dependent 
cross-section, neutron-proton cross-section is three times larger compared to neutron-neutron 
and proton-proton cross-section that will enhance binary collisions. Moreover, in case of 
neutron rich system, due to more repulsion more squeeze-out can be seen. 



To further strengthen our interpretation of the results of fig.l, we display in fig. 2, the final 
phase space of nucleons for X-Y plane and Px-Py plane for the reaction of 124 Ag4 7 + 124 Ag i7 
at an incident energy of 100 MeV/nucleon (below transition energy), 236 MeV/nucleon (at 
transition energy) and 350 MeV/nucleon (above transition energy) in upper, middle and 
below panels respectively for a randomly selected event. It has been observed that the 
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FIG. 2: Phase space distribution of nucleons in X-Y plane(left panel) and Px-Py plane (right 
panel). The reaction under study is lu Ag 4 j + 124 Ag^. The panels from top to bottom are 
representing the phase space of nucleons at different energies. 

distribution of the nucleons is ellipse like for E = 100 and 350 MeV/nucleon i.e below and 
above the transition energy and spherical for E = 236 MeV/nucleon i.e at the transition 
energy. 

To study the effect of isospin dependence of cross-section and charge asymmetry on 



dN 



,— , we display in Fig.3, the transverse momentum dependence of —7— for the reactions of 

ptapt o? r ptapt 

124 Ag 47 + 124 Ag 47 and 124 Pr 59 + 124 Pr 59 at an incident energy E = 100 MeV/nucleon. The 
figure reveal following points: 

(a) As the transverse momentum increases decrease. Which is quite obvious, because 
with increase in transverse momentum, the number of particles in that particular bin with 



7 





5.0x1 0" 3 




4.0x1 0" 3 




3.0x10 




2.0x10 


> 

CD 


2.0x1 0" 3 








1.5x10~ 3 








1.0x10" 3 


~o 




Cl 


5.0x1 0" 4 


~Z. 


3.0x10 






2.0x1 0" 3 




1.0x10" 3 




0.0 



' l ' l ' l ' l ' l ' 

— oS. FN 

ISO *v 

^noiso 


124,-. . 124,-. 

Pi" + Pr„ 


I 


i i 

E = 100 MeV/ nucleon 




I ' I i I i I i I ' I i I I 



50 100 



50 100 150 



<p t (MeV/c)> 



FIG. 3: (color online) Transverse momentum dependence of for the reactions of 124 ^4<?47 + 124 
Ag A1 (left) and 124 Pr 59 + 124 Pr 59 (right). 

large transverse momentum decreases. This shows that after the collision, momentum is 
not equally transfered among the nucleon. Some nucleon suffer hard collision while other 
suffer soft collision. 

(b) The value of is more for neutron rich system 12i Ag i7 + 124 Ag^ than neutron 

deficient system 124 Pr^g + 124 Pr$g due to increase in repulsive forces among nucleon. 



dN 
ptdpt 



is sensitive to different nucleon-nucleon cross-section. Its value is more in case of 



isospin dependent cross-section. This happens because in the case of isospin dependent 
cross-section, neutron-proton cross-section is three times large compared to neutron-neutron 
and proton-proton cross-section that will enhance binary collisions. But this increase is not 
uniform for free nucleons. 
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FIG. 4: (color online) Azimuthal angle dependence of ^cos2(j>)) ^ or ^ ree nuc l eons (upper panel) 
and LMF's (lower panel) at incident energies 50, 100, 150 and 250 MeV/nucleon. 



In fig.4, we display the azimuthal angle dependence of 



dN 



d({Cos2<j>)) 



for free nucleons and 



LMF's at incident energies 50, 100, 150 and 250 MeV/nucleon for the reaction of 1M Sn 50 + 



124 



Sn 50 . We note: 



As the energy increases, slope of the curve increases for both free nucleons and LMF's. This 
happens because, as the energy increases the thrust will also increases which will enhance the 
out-of-plane flow of the nucleon. The increase in slope is uniform in case of free nucleons but 
non-uniform in case of LMF's. Which indicates that emission of free nucleon is symmetrical 
but emission of LMF's is not symmetrical about reaction plane. This is something interesting 



which we were not expecting. 
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FIG. 5: (color online) Transverse momentum dependence of p ^ pt for free nucleons, LMF's and 
IMF's at incident energies E = 100, 150 and 250 MeV/nucleon. 

In fig. 5, we display the transverse momentum dependence of for free nucleons, LMF's 
and IMF's at incident energies E = 100, 150 and 250 MeV/nucleon for the reaction of 
124 S , ra 5 o + 124 Sn 50 . We note that as the energy increases, decreases. This happens 
because, with increase in energy the particles with large transverse momentum will decrease 
in that particular bin thus the curve shift downward for free nucleons, LMF's and IMF's. 
To further strengthen our interpretation of the results, we display in fig. 6 the energy 



dN 



for free nucleons and LMF's for the reactions of 12A Sn 5 n + 124 Sn 



dependence of d{{Cos2<t>)) 
(N/Z = 1.48) and 40 S 16 + 40 S 16 (N/Z = 1.5) for (Cos2(f)} = -1, 0, +1. One can note that, 
once the free nucleons and LMF's at (Cos2(j)) = and -1 are normalized with free nucleons 



and LMF's at (Cos2(/)) = +1 at the starting point of energy, we see that their behavior with 
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FIG. 6: (color online) The energy dependence of ^((Co^)) ^ or ^ ne reactions of 12A Sn^Q + 124 Sn^o 
and 40 S 16 + 40 5i 6 . 

respect to the energy is similar for both the free nucleons and LMF's. Although, the N/Z 
of both the reactions are nearly equal but the behavior of energy dependence of d ^os24>)) 
is different for both the reactions. One can see from the fig. 5, that out-of-plane emission 
is more in case of 124 Sri5o + 124 Sn^o than in-plane emission. But the behavior is entirely 
different for Sie + 40 5*16 where the in-plane emission is more than out-of-plane emission. 
The distribution of free nucleon is equal for in-plane and out-of-plane flow with respect to 
distribution corresponding to vanishing flow. But for LMF's the distribution is asymmetric. 
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III. SUMMARY 



Using the isospin dependent quantum molecular dynamics model, we have studied the 
effect of charge asymmetry and isospin dependent cross-section on ^tt^^h and Sim- 
ulations have been carried out for the reactions of 124 X m + 124 X m , where m = (47, 50 and 
59) and 40 S\q + 40 Siq. Our study showed that distribution of nucleons and fragments is not 
symmetric in space for in-plane and out-of-plane emission. 
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